Dirac particle in gravitational field 
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Classical dynamics of spinning zero-size objects in an external gravitational field is derived from 
the conservation law of the stress-energy and spin tensors. The resulting world line equations differ 
from those in the existing literature. In particular, the spin of the Dirac particle does not couple 
to the background curvature. As a check of consistency, the wave packet solution of the free Dirac 
equation is considered. The resulting equations are shown to include a constraint that relates the 
wave packet spin to its orbital angular momentum. In the zero-size limit, both contributions to the 
total angular momentum disappear simultaneously. 



00 

o 

O 

(N 

C 

in 

On 



I , 



> 
00 
(N 
(N 
CO 

\6 
o 
oo 
o 



X 



PACS numbers: 04.40.-b 

Introduction. The problem of particle motion in ex- 
ternal gravitational field is usually addressed by using the 
Mathisson-Papapetrou method [l], Q. One starts with 
the covariant conservation law of the stress-energy ten- 
sor of matter fields, and analyzes it under the assumption 
that matter is highly localized in space. In the lowest, 
single-pole approximation, the moving matter is viewed 
as a point particle. In the pole-dipole approximation, its 
non-zero size is taken into account. 

The results found in literature can be summarized as 
follows. Spinlcss particles in the single-pole approxima- 
tion obey the geodesic equation. In the pole-dipole ap- 
proximation, the rotational angular momentum of the 
localized matter couples to spacetime curvature, and pro- 
duces geodesic deviations [l|, 0, H, S @1 ■ If the particles 
have spin, the curvature couples to the total angular mo- 
mentum 0, SB Ed, E|- 

What we are interested in is a consistent single-pole 
analysis of spinning particles in an external gravitational 
field. This is motivated by the observation that single- 
pole approximation eliminates the influence of particle 
thickness, and allows the derivation of the pure spin- 
curvature coupling. In fact, this is the only way to see 
the influence of curvature on the spin part of the total 
angular momentum. The ambiguous algebraic decom- 
position of the total angular momentum into spin and 
orbital contributions are of no help. What we need is 
a truly zero-size object. The existing literature on the 
subject turns out not to have this sort of prediction. 

The results that we have obtained are summarized as 
follows. Trajectories of spinning zero-size massive parti- 
cles generally deviate from the geodesic lines. The devi- 
ation is due to the spin-curvature and spin-torsion cou- 
plings. These turn out to be different from what has been 
believed so far. In particular, the spin of the free Dirac 
point particle does not couple to the curvature. If it is 
viewed as a wave packet solution of the Dirac equation, 
it does not couple to the torsion either. The Dirac point 
particles behave as spinless objects. 

In what follows, we shall outline our derivation, and 
make the necessary comments. We begin with the covari- 
ant conservation of the fundamental matter currents 



stress-energy tensor t p v and spin tensor er A M „: 
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ant derivative with the nonsymmetric connection P M „ 
The second equation clearly shows that the antisymmet- 
ric part of the stress-energy tensor is not a dynamic vari- 
able. Indeed, given the spin tensor a Xp ' lJ , the antisymmet- 
ric part is completely determined. We can use the 
equation (jlb[) to eliminate components, and rewrite 
the conservation equations in terms of -r^' and a x>iv 
alone. 

Let us briefly explain how the assumption of highly 
localized matter simplifies the conservation equations (fT]) . 
First, we expand the stress-energy and spin tensors of the 
localized matter in multipoles around a suitably chosen 
timelike line x^ = z p (t). Using the results of Ref. 
a spacetime function F(x) is written as a series of 8- 
function derivatives: 

8^{x 



F(x) = / dr 
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The line parameter r is chosen to be the proper distance, 
dr 2 = g^ v dz p dz u , and V p = D p (r -> {}) is the Rie- 
mannian covariant derivative. Here, ( A ) stands for the 
Levi-Civita connection. The expansion (|2|) is manifestly 
covariant, with the multipole coefficients M Pl "- Pn being 
tensors with respect to spacetime diffeomorphisms. The 
adopted metric signature is +2, and g = det^g^) is the 
metric determinant. 

The decomposition is suitable for treating matter 



which is well localized around the line 



z^(t). In 



fact, the variable F(x) must drop exponentially to zero 
as we move away from the line if we want the series ([2]) 
to be well defined. If this is the case, the coefficients 
M pi --- Pn get smaller as n gets larger. In the lowest, single- 
pole approximation, all higher M's are neglected. In the 



2 



pole-dipole approximation one neglects all but the first 
two M-coefficients. 

World line equations. We expand the stress-energy 
and spin tensors of the localized matter in multipolcs 



around the line 



'(r). Being concerned with the 



case of point particles, we choose to work in the single- 
pole approximation. Thus, 



Xfiu 



„, ,6W (x-z) 



(3a) 



(3b) 



This definition of single-pole approximation differs from 
what is found in the existing literature [H, [l(J. There, 
the antisymmetric part of the stress-energy tensor 
had an equal treatment as t^ v '. As is not an inde- 
pendent variable, this imposed unnecessary constraints 
on fr AAIly . In particular, the spin of the Dirac particle 
was ruled out. To overcome this problem, the authors of 
Ref. [l(| abandoned single-pole in favour of pole-dipole 
approximation. Their subsequent limit of vanishing or- 
bital angular momentum brought them a non-trivial spin- 
curvature coupling, even in the Dirac case. In what fol- 
lows, we shall demonstrate that it is not quite so. 

Let us substitute ([3]) in Jl}, and solve the obtained 
equations in terms of the unknown variables z M (r), c pv (r) 
and c Xpv {r). The procedure is very much the same as 
described in Ref. [121 ]. First, we eliminate the equation 
(llb[) , and rewrite (flaj) in terms of the of the independent 
variables t^"' and a Xp - v . We then multiply (fl"aj) with 
an arbitrary vector valued function f^ix), and integrate 
over the spacetime. The resulting equation depends on 
the function / M and its first and second covariant deriva- 
tives, evaluated on the line x^ — z^(t). Owing to the ar- 
bitrariness of the function fp,(x), the terms proportional 
to its independent derivatives separately vanish. The ob- 
tained equations are manifestly covariant with respect to 
spacetime diffeomorphisms. 

Skipping the details of the calculation, which will be 
published elsewhere, we present here the final result. It 
consists of two sets of world line equations. The first is 
a set of algebraic equations which determine the coeffi- 
cients cP" and c Xpv in terms of the free parameters m(r), 
s^(t) and s A ^(r): 



(4a) 
(4b) 



Here, u p = dz^/dr is the particle 4-velocity, K^xv = 
(%fj,\ — TfiXv — T Xvp )/2 is the contorsion tensor, and V 
stands for the Riemannian covariant derivative along the 
particle trajectory (Vi^ = dv p /dr + {£ p }v x u p ). The pa- 
rameters s pv and s Xpv are totally antisymmetric tensors: 
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The algebraic equations (|4]) differ from what we find in 
literature by the form of (|4b[) . 

The second set of equations are the differential equa- 
tions that determine the particle world line and spin pre- 
cession. They read: 



Vp" - u»s<> x R» vpX 
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where 



r^vpX 



(5a) 
(5b) 



p p = mw 11 + 2u p (Vs' ip + D MP ) , 
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-K Xp ^c^ px . 



Here, the symbol R^ up \ stands for the Riemann curva- 
ture, R^vpX = Tl^vpxir — > {}), and both external fields, 
R^vpX and K pXv , are evaluated on the line x p = z p (t). 

The obtained single-pole equations differ from the 
known pole-dipole result 0, H, @, 2, U | by the presence 
of the constraint (|4b[) . It is a consequence of our assump- 
tion that the particle has no thickness, and therefore no 
orbital degrees of freedom. In the existing literature, an 
analogous but more restrictive constraint appears in this 
regime [H, [l(| • This is because the antisymmetric part of 
the stress-energy tensor has been treated as an inde- 
pendent variable, in spite of the restriction (jlbj) . In our 
approach, the only independent variables are a Xpv and 
t^"' , and the resulting constraint (|4bj) is not so strong. 
In particular, it does not rule out the free Dirac field, or 

any other massive elementary field. Indeed, the formula 
c Xpv _ u x s pv + \ u \p, s v\x for tlle gpin tensor f t ^ e e j e _ 

mentary particle of spin s (see Ref. [HI]) is a special case 
of (TO- 

In what follows, we shall examine the special case of 
spin 1/2 pointlike matter. Surprisingly, we shall discover 
that spin 1/2 docs not couple to the curvature, leading 
to geodesic trajectories in torsionless spacetimes. 

Dirac particle. The basic property of Dirac matter 
is the total antisymmetry of its spin tensor a Xpv . As a 
consequence, the coefficients c Xpv are also totally anti- 
symmetric, and the constraint (|4b() implies 



0. 



The vanishing of the s Ml/ component of the spin tensor 
has far reaching consequences. First, we see that the 
spin-curvature coupling disappears from the world line 
equation (|5a[) . Second, the spin precession equation (|5b[) 
becomes an algebraic equation. If we define the spin 
vector s p by s pvp = e pl/pX s x , and the axial component 
of the contorsion K p as K p = e^ pX K upXl the equations 
© become 



V (m/ + K^s^u^ + ^s v VK v = , (6a) 
. (6b) 



K [ lsi 
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Here, e^px = \/—g£pvpX is the covariant Levi-Civita 
symbol, and the decomposition of a vector into compo- 
nents parallel and orthogonal to the world line is used: 
= Vj_ + uit M , v^_Up = 0. As we can see, the spin 
couples only to the axial component of the contorsion, 
which means that Dirac point particles follow geodesic 
trajectories in torsionless spacetimes. At the same time, 
the absence of torsion trivializes the equation (|6b[) . and 
no information on the behavior of the spin vector is avail- 
able. 

To make it clear, only zero-size Dirac particles lack 
the spin-curvature coupling. In the case of thick par- 
ticles (pole-dipole approximation), curvature couples to 
the total angular momentum, thus to both, spin and or- 
bital angular momentum. Therefore, our result does not 
compromise the known pole-dipole result. It just explains 
its behaviour in the zero-size limit. 

If the background torsion has nontrivial axial compo- 
nent, a geodesic deviation appears, but also a very strong 
constraint on the spin vector. Indeed, the equation (|6bj) 
implies that the orthogonal component of s M always ori- 
ents itself along the background direction K± . This un- 
usual behavior suggests that the spin vector of the Dirac 
point particle might be zero, after all. In fact, the world 
line equations ^ and ([5]) are derived under very gen- 
eral assumptions of the existence of pointlike solutions in 
an arbitrary field theory. They do not care about pecu- 
liarities of specific theories, or specific types of localized 
solutions. What might happen in a particular theory 
is that additional constraints appear to ensure the ex- 
istence of localized solutions. In what follows, we shall 
analyze the wave packet solutions of the flat space Dirac 
equation, with the idea to check if they can be viewed as 
point particles. We shall see that this is possible only if 
the particle spin goes to zero. 

Wave packet. Wave packet is a solution which is well 
localized in space, but resembles a plane wave inside. To 
be viewed as a particle, its size I is considered in the 
limit £ — * 0. At the same time, the particle stability is 
achieved in the limit X/£ — > 0, where A is its wavelength. 
We shall consider the free Dirac equation, and construct 
wave packets as follows. In the initial moment t — 0, we 
choose the configuration ^(r,0) = g(r)ijj p (r, 0), where 
ip p (x) = a(k) exp(ifc M x M ) is a plane wave solution of the 
Dirac equation, and g(r) ~ exp(— r 2 /£ 2 ) is a function 
that cuts out a small piece of the plane wave. Using 
the Dirac equation, we can calculate time derivatives, 
and thereby determine its time evolution. In fact, we 
only need first time derivatives, as neither t* 1 " nor a x ^ v 
depend on higher derivatives: 

t> = i^p d»ip , o- x ^ = £ A ^^ 7 57p^ • (?) 

In addition, there is a conserved current j^ = tpj 11 ^, 
associated with the extra U(l) symmetry of the Dirac 
action. The wave packet expressions of these currents 



are obtained straightforwardly, but we shall skip these 
details. Instead, we shall discuss two interesting relations 
that we have discovered to further constrain the currents. 
The first, 

VmO' A = 0, (8) 

shows that the spin vector is orthogonal to u^, as 
expected. This follows from the fact that j 71 oc w' 1 in the 
single-pole approximation. The second relation, 

r^jVocj", (9) 

has even more striking consequences. Its physical mean- 
ing is best seen in the rest frame where it takes the form 
r a0 = j a = 0. It tells us that the two currents are mutu- 
ally proportional, i.e. that the energy and charge flow in 
the same direction. In the limit X/£ — > 0, one can show 
that it implies the constraint 

where div stands for the total 3-divergence. In the single- 
pole approximation, only the space integral of this rela- 
tion is needed. Then, the divergence term disappears, 
and we end up with the constraint that relates the wave 
packet spin to its orbital angular momentum. As orbital 
degrees of freedom are neglected in this approximation, 
so is the spin. Thus, the particle orbital angular momen- 
tum and spin disappear simultaneously in the limit £ — > 0, 
X/£ -f 0. 

To summarize, we see that the spin vector s M vanishes 
in the single-pole approximation, and the particle tra- 
jectory becomes a geodesic line even in the presence of 
torsion. This is a valid conclusion only if Dirac point par- 
ticles are viewed as zero-size wave packets. It does not ap- 
ply to other localized field configurations. For example, 
one could consider a localized bound state of Dirac and 
electromagnetic fields. Its spin-tensor is still totally an- 
tisymmetric (which implies the vanishing spin-curvature 
coupling), but nothing proves its disappearance in the 
zero-size limit. The validity of the above conclusion also 
requires some sort of equivalence principle to hold. This 
is because the considered wave packet is a solution of 
the free Dirac equation, and the inclusion of curvature 
or torsion may destroy it. What we can do is to con- 
sider weak gravity, so that terms quadratic in curvature 
and torsion are neglected. In that case, the free wave 
packets are a good approximation to the exact solution, 
which implies that Dirac point particles behave as spin- 
less objects in a weak gravitational field. They can still 
probe the spacetime curvature, but for the probe of the 
background torsion, one needs a thick particle. 

Let us note here that one could use the alternative 
WKB approach to study Dirac field in curved back- 
grounds. In such an approach, the starting point is the 
covariant Dirac equation rather than the general conser- 
vation laws. However, whatever solution of the covari- 
ant Dirac equation is used, it must obey the covariant 



conservation equations (JTJ) . This is because the conser- 
vation laws |T]) are a general property of all covariant 
Lagrangians, the covariant Dirac Lagrangian being just 
one of those. Once we check that the same approximation 
is used in both approaches, the two predictions must not 
contradict each other. In particular, if the wave packet 
solution of the covariant Dirac equation is considered in 
the zero-size limit, its trajectory will belong to the class 
of trajectories defined by our world line equations. The 
WKB analysis of the covariant Dirac equation has al- 
ready been considered in literature (see, for example, Ref. 
[HI). However, this kind of analysis is not suitable for 
the zero-size limit, as the considered wave configuration 
is far from resembling a point-particle. What is called 
particle trajectory in this approach is rather a ray of ge- 
ometric optics. 

Let us close our exposition with the speculation that 
Dirac particle is not the only one that possesses these 
features. Indeed, the relations analogous to (J9j) might 
exist quite generally. There is nothing special about the 
statement that all the particle charges flow in the same 
direction. This is something one would expect to hold for 
any type of localized matter. However, the consequence 
(fT0| has been derived for the Dirac wave packet solution 
only. If this turns out to be the general property of wave 
packet solutions, the disappearance of spin in the zero- 
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size limit might be a feature of all massive point particles. 
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